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We study random fiber bundle model (RFBM) with different threshold strength distributions 
and load sharing rules. A mixed RFBM within global load sharing scheme is introduced which 
consists of weak and strong fibers with uniform distribution of threshold strength of fibers having 
a discontinuity. The dependence of the critical stress of the above model on the measure of the 
discontinuity of the distribution is extensively studied. A similar RFBM with two types of fibers 
belonging to two different Weibull distribution of threshold strength is also studied. The variation 
of the critical stress of a one dimensional RFBM with the number of fibers is obtained for strictly 
uniform distribution and local load sharing using an exact method which assumes one-sided load 
transfer. The critical behaviour of RFBM with fibers placed on a random graph having co-ordination 
number 3 is investigated numerically for uniformly distributed threshold strength of fibers subjected 
to local load sharing rule, and mean field critical behaviour is established. 
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I. INTRODUCTION 

Sudden catastrophic failure of structures due to unex- 
pected fracture of component materials is a concern and 
a challenging problem of physics as well as engineering. 
The dynamics of the failure of materials show interesting 
properties and hence there has been an enormous amount 
of study on the breakdown phenomena till nowi^^i^**^. 
A total analysis of the fracture of materials should in- 
clude the investigation of a possible fracture before the 
occurrence of large deformations and hence the determi- 
nation of the critical load beyond which if the load on 
the system increases, the complete failure occurs. Future 
progress in the understanding of fracture and the applica- 
tion of that knowledge demands the integration of contin- 
uum mechanics with the scientific disciplines of material 
science, physics, mathematics and chemistry. The main 
challenge here is to deal with the inherent stochasticity 
of the system and the statistical evolution of microscopic 
crack to accurately predict the point of final rupture. 

Having said about the complexity involved in fracture 
processes, one can still get an idea of fracture phenomena 
by studying simplified models. The simplest available 
model is Fiber Bundle Model(FBM)££&. It consists of N 
parallel fibers. All of them have their ends connected to 
horizontal rods as shown in Fig. 1. The disorder of a real 
system is introduced in the fiber bundle in the form of 
random distribution of strength of each fiber taken from 
a probability distribution P(cr) and hence called Random 
Fiber Bundle Model (RFBM). The strength of each fiber 
is called its threshold strength. As a force F is applied 
externally on a bundle of N fibers, a stress a = F/N acts 
on each of them. The fibers which have their thresh- 
old strength smaller than the stress generated, will break 
immediately. The next question that arises is the affect 
of breaking of fibers on the remaining intact fibers i.e., 
one has to decide a load sharing rule. The two extreme 
cases of load sharing mechanisms are Global Load Shar- 
ing (GLS)£2 and Local Load Sharing (LLS) 10 i n i 12 . In 
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FIG. 1: The vertical parallel lines correspond to fibers which 
are connected to two horizontal plates (shown with bold 
lines). The upper plate is fixed and the lower plate is mov- 
able. Fig. 1(a) shows a smaller force Fi applied on the bundle 
producing a strain S. Fig. 1(b) corresponds to a higher value 
of load, leading to a strain S and fibers begin to break. 

GLS, the stress of the broken fiber is equally distributed 
to the remaining intact fibers. This rule neglects local 
fluctuations in stress and therefore is effectively a mean 
field model with long range interactions among the ele- 
ments of the systemic. On the other hand, in LLS the 
stress of the broken fiber is given only to its nearest sur- 
viving neighbours. It is obvious that the actual breaking 
process involves a sharing rule which is in between GLS 
and LLS. Several studies have been made which considers 
a rule interpolating between GLS and LL S 14 ' 1 ^ . 

For a given force F, the fibers break and distribute their 
load following a load sharing rule until all the fibers have 
their threshold strength greater than the redistributed 
stress acting on them. This corresponds to the fixed point 
of the dynamics of the system. As the applied force is 
increased on the system, more and more fibers break and 
the additional stress is redistributed amongst the other 
fibers. There exist a critical load (or stress a c ) beyond 
which if the load is applied, complete failure of the sys- 
tem takes place. Most of the studies on FBM involve the 
determination of the critical stress a c and the investiga- 
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tion of the type of phase transition from a state of partial 
failure to a state of complete failure. It has been shown 
that a bundle following GLS has a finite value of critical 
stress and belongs to a universality class with a specific 
set of critical exponents^^ whereas there is no finite crit- 
ical stress <7 C at thermodynamic limit in the case of LLS 
in one dimensioniflii£iiSii£. On the other hand, LLS on 
complex network has been shown to belong to the same 
universality class as that of GLS with the same critical 
exponents^. 

The paper is organised as follows. Section II consists 
of mixed fibers with uniform distribution (UD) of thresh- 
old strength of fibers and GLS. The critical behaviour of 
such a model is explored analytically. Mixed fibers with 
Weibull distribution (WD) of threshold strength along 
with GLS is studied using a quasistatic approach of load 
increase in section III. One dimensional, one sided LLS 
with strictly uniform distribution (UD) is presented in 
section IV. Section V consists of a bundle of fibers whose 
tips are placed on a random graph with co-ordination 
number 3. The main conclusions are presented in section 
VI. 



II. MIXED FIBERS WITH UNIFORM 
DISTRIBUTION AND GLS 



Let us consider a 'mixed' fiber bundle model with 
normalised density distribution of threshold strength as 
shown below: 
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FIG. 2: Mixed Uniform Distribution. 
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The present model clearly consists of fibers of two differ- 
ent types. This classification is based on whether a ran- 
domly chosen fiber from the bundle belongs to class A of 
fibers (0 < crth < ci) or class B of fibers {a 2 < crth < 1) 



as shown in Fig. 2. In other words, the model consists 
of weak and strong fibers but no fibers with strengths 
between ax and 02, which is the forbidden region for the 
fiber strength. Since we are considering GLS here, the 
exact location of class A or class B fibers is immaterial. 
The aim is to determine the critical stress cr c for this 
model such that if the external stress exceeds cr c , the 
bundle breaks down completely. 

If a fraction x of the total fibers belong to class A and 
the remaining 1 — x to class B , then to maintain the 
uniformity of the distribution(Eq. (1)), we must have 
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The above equation provides the relationship between x, 
ax and a 2 showing that once x and ax are set, 02 is fixed 
automatically. This equation at the same time puts a 
restriction on ax- One can check that if ax > x, 02 
becomes smaller than ax, which is impossible with the 
given form of distribution (Eq. (1)). We shall therefore 
restrict the value of ax to be smaller than x. 

The mechanism of the failure of RFBM is the follow- 
ing: under an external stress cr, a fraction of fibers having 
threshold strength less than the applied stress fail imme- 
diately. The additional load due to breaking of these 
fibers is redistributed globally among the remaining in- 
tact fibers. This redistribution causes further failures. 
Let us define Ut = Nt/N as the fraction of unbroken 
fibers after a time step t. Then, 
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The recurrence relation between Ut and Ut+x for a 
given applied stress a is obtained as^: 



Ut+x = 1 - P(tr t ) = 1 



(3) 



This dynamics propagate (in discrete time) until no fur- 
ther breaking takes place. 

If the initial applied load is very small so that the re- 
distributed stress is always less than 02, then no fibers 
from class B will fail. Thus for breaking of class B fibers, 
the redistributed stress has to be greater than a 2 - Let 
at (the redistributed stress) cross a% at some instant t. 
Then the critical stress of the system can be calculated 
as follows: 
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At fixed point of the dynamics, when no further break- 
ing takes place, Ut+x — U t = U* and we get, 
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U* = 1 - [ 
giving, 
U* 
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It can be shown easily that the solution with + in the 
parenthesis gives a stable solution. Therefore, a c is the 
critical stress of the system because if the applied stress 
a is less than a c , the system reaches a fixed point. If the 
applied stress is greater than a c , U* becomes imaginary 
and the bundle breaks down completely. It is interesting 
to note that at a — a c , U* = 1/2[1 — ((72 — eri)] i.e., more 
than half of the fibers are still unbroken at critical point. 

We should now look at the restrictions on a\, a% and 
x (other than <j\ < x) so that the above calculations are 
valid. From Eq. (4), we find that if the width of the for- 
bidden region (a 2 — <Ji) is greater than half, U*(a c ) > 1 
which is unphysical. Therefore, (a 2 — a\) should be 
smaller than 1/2 for all values of x for the above cal- 
culations to hold. Using this restriction in Eq. (5), we 
find that the critical stress a c will always be less than 
1/2. Now, looking at the physical meaning of U* , we 
must also check that the fraction of broken fibers at the 
fixed point corresponding to critical stress a c must be 
greater than the fraction of fibers having their threshold 
stress smaller than the critical stress. We summarise all 
the restrictions below: 

01 <x, (7 2 - en < 0.5 and (1 - U* (ct c )) > P(cr c ) (6) 

In the table below, we enlist some allowed values of 
(7i, (72 for different x and also the corresponding critical 
stress. 
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From the table, we observe that the critical stress de- 
pends entirely on the width of the forbidden region, i.e., 
<7 2 — 0"i • Also, when o\ is increased for a fixed x, 02 de- 
creases and hence a c decreases. Physically, decreasing a 2 
(with x fixed) means that lowest threshold of the strong 



fibers (class B) gets lowered. It is found that a\ cannot 
be very small compared to x, else all the above mentioned 
conditions do not get satisfied. 

One can define an order parameter as shown below: 



2[1 - ((72 - <7i)]ET - 1 = (a c - a)* = (a c - a) 
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The order parameter goes to as a — > a c following a 
power law (a c — <r)' . Susceptibility can be defined as 
the increment in the number of broken fibers for an in- 
finitesimal increase of load. Therefore, 

dm 

X = -r- where m = N[l-U*(a)] 
da 

Hence, 

X oc (cr c - a)~ 



(a c - a)' 



(8) 




FIG. 3: Variation of a c with <Ji - 
a distribution as in Fig. 2 



<7i for a mixed RFBM with 



As expected, the model follows mean field exponents 
with (3 = 1/2 and 7 = 1/2. The interesting features of 
this model are: In this model, the forbidden zone has a 
nontrivial role to play. (i)By changing (02 — 01) and x, we 
can tune the strength of the system. However there are 
restrictions on the parameters as discussed above. a c as a 
function of ((72 — fi) is shown in Fig. 3. (ii) If x = 0, (72 7^ 
0, (7i = 0( from Eq. 2) which gives a c = 1/4(1 — (72) and 
we observe an elastic to plastic deformation as discussed 
in reference 9. In this case, for a < 02, the fiber in 
the bundle obey Hook's law. As soon as a > 02, the 
dynamics of failure of fibers start. (iii)If (72 — (7i=0, the 
problem is that of a simple UD with threshold strength 
between and 1 which has a critical stress a c = 0.2E 9 . 

These results confirm that the failure of mixed fiber 
bundle model with Global Load Sharing belongs to a 
universality class with critical exponents (3=1/2 and 
7=1/2- 



III. MIXED FIBERS WITH WEIBULL 
DISTRIBUTION AND GLS 

In the Weibull distribution of threshold strength of 
fibers, the probability of failure of each element has the 
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form 
P(a) = 1 



-(—)" 



(9) 



where <tq is a reference strength and p is the Weibull 
index. We now consider a mixed RFBM as in section II 
with WD of threshold strength of fibers. As in the UD 
case (of Sec. (2)), here also a fraction x of fibers belong 
to the class A (WD with p = 2) and the remaining (1 — x) 
fraction belong to the class B (WD with p = 3) with the 
reference strength <tq equal to 1 in either case. 

For the conventional WD, consider a situation where 
the stress on each fiber increases from g\ to a 2 . The 
probability that a fiber randomly chosen from the WD 
survives from the load o\ but fails when the load is 02, 
is given by 



p(ai,<r 2 ) = 



= i- e -K-<). 



Thus the probability that the chosen fiber that has sur- 
vived the load ai also survives the (higher) load a% is 



9(0-1, er 2 ) = e 



The key point is that the force 



F on the bundle is increased quasistatically so that only 
one fiber amongst the remaining intact fibers break. As 
mentioned previously, the dynamics of breaking will con- 
tinue till the system reaches a fixed point. The process 
of slow increase of external load is carried on up to the 
critical stress er c . 

To deal with the nature of the present 'mixed' random 
fiber bundle model, the method used by Moreno, Gomez 
and Pacheco^i is generalised in the following manner. It 
is important here to keep track of number of unbroken 
fibers in each distribution separately. Let us assume that 
after a loading is done and a fixed point is reached, N k2 
and N k3 are the number of unbroken fibers corresponding 
to p = 2 and p = 3 distribution respectively, where each 
fiber has a stress a k . 
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FIG. 4: Variation of fraction of broken fibers with external 
load per fiber a for a RFBM with two distributions corre- 
sponding to p = 2 and p = 3. N=50000, :r=0.5 

We define q 2 and as follows: 



03(0-1, 02 J = e ^ 

The load (N k2 + N k3 )ai that has to be applied to break 
one fiber is given by <j\ — min(ai 2 ,ai 3 ) where ct/ 2 (or 
cr/ 3 ) is the next weakest fiber in the p = 2 (or p = 3) 
distribution of strength of fibers and is obtained by the 
solution of the following equations. 

N k2 ~ 1 = N k2 q 2 (cr k ,ai 2 ) 
and N k3 - 1 = N k3 q 3 (a k , ai 3 ) 

which gives 
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The breaking of one fiber and the redistribution of its 
stress causes some more failures. Let us say that dur- 
ing this avalanche, at some point before the fixed point 
is reached, there are N k2 and N ka number of unbroken 
fibers belonging to the two distributions where each fiber 
is under a stress a%. This stress causes some more fail- 
ures and as a result N k2 and N k3 fibers are unbroken. 

Let the new stress developed is a k . The number of fibers 
which survive 5> and a k are obtained using the relation 



N k 2 = N k2 q 2 (a k ,a k ) 



(12) 
(13) 
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The stress on each fiber is now equal to (N k + 

<K/(< + <)• E q (12) and Eq. (13) are used 
again and again until a fixed point is reached. The fixed 
point condition is given by (N k2 + N k3 ) — {N k2 + N k3 ) < e 
where e is a small number (0.001). However, the critical 
behaviour does not depend on the choice of e. 

After the fixed point is reached, stress 07 is calculated 
once again as mentioned before and the whole process is 
repeated till complete failure occurs at a c . The above 
mentioned approach has the advantage that it does not 
require the random averaging of Monte Carlo simula- 
tions. 

Fig. 4 shows the fraction of total number of broken 
fibers as a function of applied stress a for the case 
x = 0.5. The graph clearly shows the existence of a 
critical stress, er c = 0.46, which lies between the criti- 
cal stress of p = 2 (0.42) and p = 3 (0.49) distribution 
. These theoretical values are obtained from the analyt- 
ical expressiorAjc = (pe)^ 1 ^. The avalanche size S is 
defined as the number of broken fibers between two suc- 
cessive loadings. It diverges near the critical point with 
an exponent 7 = 1/2 as (a c — o)~ 7 - Scaling behaviour 
of avalanche size S is shown in Fig. 5. This confirms 
the mean field nature of GLS models even with mixed 
fibers. The most attractive feature of the mixed RFBM 
is shown in Fig. (6) which contains the variation of a c 
with the fraction x of class A fibers. Interestingly, a c 
decreases linearly as x is increased. 
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FIG. 5: Scaling behaviour of avalanche size as the critical 
point is reached for the mixed model. Also shown is a straight 
line (dotted) with slope (-1/2). Here, N=50000 and x = 0.5. 
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FIG. 6: Variation of critical stress a c with x. Slope of the 
straight line is 0.068. N=50000 



IV. ONE DIMENSIONAL, ONE SIDED LLS 

In numerical studies, it is difficult to generate a strictly 
uniform distribution of the threshold strength of random 
fibers. In this work, we have addressed the problem of 
random fibers with UD of threshold with LLS scheme 
using an exact treatment of one sided load distribution 10 . 
The complication arising in the present approach due to 
the uniform nature of the threshold distribution will also 
be emphasised later. 

In one sided LLS, the additional load due to breaking 
of a fiber is transferred only to its neighbour in a given di- 
rection. The fibers are placed along the circumference of 
a circle with periodic boundary conditions. The method 
is based on the k- failure concept. If k is chosen appropri- 
ately, then more than k adjacent failures is equivalent to 
complete failure of the system. The scheme involved is 
as follows: let pi = p{ia) is the probability of failure of a 
fiber supporting a load ia (where a is the initial applied 
load per fiber). This is same as the probability of failure 



of all the previous fibers. Let Pjv(c) represents the 

probability that a system of N fibers break when an ini- 
tial force equal to Na is applied and Qn{&) = 1 — P/v(c) 
is the probability that the bundle does not break because 
of initial load Na. We have defined the critical load a c 
by setting Pm{u c ) — 0.5. However, any other choice of 
Pn(& c ) does not alter the result. Qn{&) can be easily 
identified as the probability that out of N fibers, there is 
no chain of adjacent broken fibers of length greater than 
k. 

Qn{v) can be calculated by using a simple iterative 
procedure described in Fig. 7. The figure shows all possi- 
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FIG. 7: Possible breaking configurations of a set with N=4 
elements when the maximum allowed crack length is k=2 for 
one sided one dimensional LLS 

ble states of a fiber in a system of 4 fibers with maximum 
number of allowed adjacent broken fibers equal to 2 (i.e., 
k = 2). The notation means the following: the symbols 1 
and 2 mean that in the corresponding position, the fiber 
breaks with a stress la, 2a respectively. l',2',3' means 
that the corresponding element does not break under the 
stress lcr, 2a and 3er respectively. This simple notation 
can be extended to a system of N fibers. If a fiber site 
has a symbol to, the two possibilities for the next fiber 
are (m + 1) and (m + 1)'. If the site has to', the two pos- 
sibilities are 1 and 1' for the next site as to' means that 
the fiber has not broken and hence no additional stress 
is given to the next fiber. But if m = k, then the only 
allowed possibility for the fiber at the next site is (k + 1)'. 
This is to ensure that we are looking at only those config- 
urations in which the complete breakdown does not take 
place. The aim is to write Qat(ct) in terms of Qjv-i(c)- 
Let Bi(N), where i < k is the sum of probability contri- 
butions of the paths of length N ending in the symbol i. 
Also, Bi'(N) with i' < k+ 1, is the sum of contributions 
of all paths ending in symbol i'. The following relations 
can then be obtained 

fe+i 

B 1 (N)=p 1 J2B i >(N-l) (14) 
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B i (N)= Pi B l _ 1 (N -1), l<i<k 

k+l 

B 1 ,(N) = q 1 ^B i ,{N-l) 

i=l 

B i ,(N)=q i B i -x(N-l), l<i <k+l. 
Thus Qn can be defined as 

k k+l 
Q N = J2Bi(N) + J2 B i'W (15) 

i=l i=l 

and hence the critical point can be obtained by demand- 
ing that Qjv(c) is 0.5 as the critical point is reached. 
B(N) can be written in a column matrix form where 
Bi(N) and Bi'(N) are its (2k+l) components. The set 
of equations represented by Eq.(14) can then be written 
as: 

B(N) = M.B(N - 1) (16) 

M can be identified as a square matrix with (2k + l) 2 
components, all independent of N. M consists of the fol- 
lowing elements: 

first row M^i = 0, M 1 ; / = pi 

second row i = 0, M 1 ' = qi 
rest of the rows: 

M hi = Sj-i^pj, M jt , = 0; 

(i,j>l) .W,- , Sj .1/,,. I) 

where i=l,...k and i =1,2 , (k + l) . 

Eq. (16) can now be written as: 

B(N) = M N B(l) (17) 

where vector B(l) is given by Bi(l) = 5i ipi, and 

5/(l) = <V,i'9i 

Since the model is based on the fc-failure concept, the 
next question that arises is what an appropriate value 
for k should be (which fixes the matrix size). This can 
be obtained by a simple criterion of convergence: if re- 
sults do not vary when matrix size is increased from k 
to k+l, then that value of k is sufficient for /c-failure 
concept. This is equivalent to saying that the probabil- 
ity of breaking when the stress on a fiber is (k + 1 )a is 1 
(by the definition of fc-failure)and hence there is no fur- 
ther change as one increases A: (as the complete failure has 
already occurred). 

The above described convergence can be clearly seen 
if one uses WD of threshold strength' 10 . As mentioned 
before, the situation is different in the case of UD. The 
method used above allows the probability of failure to 
go beyond 1 if UD is used. The probability that a fiber 
breaks when the stress is xa is xa, where x should vary 




FIG. 8: The variation of a c with k for UD is shown in Fig 8(a). 
Fig. 8(b) shows the variation of l/a c with log2N. 

from 1 to k. It may happen that for some x and a, the 
matrix elements become greater than 1 which is clearly 
unphysical as each matrix element corresponds to a prob- 
ability. It is observed that the value of a c increases with 
k and after attaining an almost constant a c , it starts 
decreasing. This decrease actually corresponds to the 
situation when the matrix elements at some stage of the 
calculation is either greater than 1 or less than and 
hence not a saturation but a peak is obtained. Fig. 8(a) 
shows the variation of a c with k where a peak is observed. 
We have taken k = 8 for our further calculations. 

In Fig. 8(b), we have shown the variation of l/a c 
with log2-/V. Clearly it is a straight line which matches 
with the previously obtained analytical and numerical 
resulta 10 i 17 i 18 . The slope of the line is 0.81. The critical 
stress vanishes in the thermodynamic limit indicating the 
absence of a critical behaviour in one sided, one dimen- 
sional RFBM. 



V. FIBERS ON A RANDOM GRAPH 

Construct a random graph having N sites such that 
each site has exactly 3 neighbours. The algorithm used 
is as follows: label the sites by integers from 1 to N where 
N is even. Connect site itoi + 1 for all i such that site N 
is connected to site 1. We call this connection as direct 
connection. Thus we have a ring of N sites. Now connect 
each site randomly to a unique site other than i + 1 and 
i — 1, thus forming N/2 pairs of sites. This connection is 
called random connection. Each site in such a graph has 
a coordination number 3. In this construction, all sites 
are on the same footing (See Fig. 9). The tip of the fibers 
are placed on these N sites and are associated with a UD 
of threshold strength. 

The system is then exposed to a total load Na such 
that each vertex carries a stress a(i) = a. If ath{i) < 
at any site i, the fiber breaks and its stress is equally dis- 
tributed to its surviving neighbours. The above condition 
is examined for each site and is repeated until there is no 
more breaking. Increase the total load Na to N(a + da). 
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FIG. 9: An example of a random graph with coordination 
number 3 and 16 sites 
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FIG. 10: The variation of ir with a for different values of 
N. Here, iV = 2 13 , 2 14 , 2 15 , 2 16 , 2 17 , 10 5 , 2 * 10 5 . The lower a c 
corresponds to a higher value of N. 



The new total load is then shared by the remaining un- 
broken fibers and once again the breaking condition is 
checked at each site. To make sure that the total load 
is conserved, if a fiber encounters a situation in which 
one of its neighbours is broken, it gives its stress to the 
next nearest surviving neighbour in that direction. This 
rule however is not applicable to the randomly connected 
fiber i.e., if the randomly connected neighbour of a bro- 
ken fiber is already broken, the stress then is given only to 
its directly connected neighbours. Once again, the load is 
increased and the process is repeated till complete break- 
down of the system. In the numerical simulations, da is 
taken to be 0.0001 and the average is taken over 10 4 to 
10 5 configurations. 

Fig. (10) shows the variation of tt, the fraction of un- 
broken fibers, with stress a for different values of N. 
Clearly, there exist a non-zero value of critical stress even 
for large N. 

To determine a c precisely, we use the standard finite 



size scaling method^ with the scaling assumption 

n(a,N) =N- a f((a-a c )N$) (18) 

where f(x) is the scaling function and the exponent v 
describes the divergence of the correlation length £ near 
the critical point, i.e., £ oc \<j — o c \~ v . If 7r oc (a c — a) 13 , we 
obtain the relation a — (3 jv. The plots of 7riV a with a for 
different N crosses through a unique point a = <r c =0.173 
with the exponent a = 0.5 (See Fig. 11). 
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FIG. 11: Precise determination of <r c from the finite size scal- 
ing form (18). The <r c at which the curves intersect is 0.173 

With this value of a c and a, one can again use the scal- 
ing relation (18) to determine the exponent v by making 
the data points to collapse to an almost single smooth 
curve as displayed in Fig. (12). The value of v obtained 
in this manner is equal to 1 which automatically gives the 
value of the other exponent f3— 0.5. This matches with 
the mean field exponents derived analytically^. A ran- 
dom graph with connections as described above is equiva- 
lent to a Bethe lattice with co-ordination number z (= 3 
in present case) for large N24. It is observed that in 
this model, accurate value of cr c (Fig. 11) and data col- 
lapse(Fig. 12) obtained by finite size scaling method is 
applicable in the limit of large N. It may be because of 
the resemblance of the underlying structure of the sys- 
tem with the Bethe lattice that it has mean field expo- 
nents. This should also be compared to Ising models with 
infinite-range interaction (which corresponds to GLS in 
the present context) 1 ^ and with short-range interaction 
on a Bethe lattice both showing mean field critical be- 
haviour. 



VI. CONCLUSION 

We have introduced a mixed RFBM with UD in one 
case and WD in the other within GLS scheme. The mixed 
RFBM with UD has a forbidden region 02 — ci so that 
the distribution of the threshold strength is discontin- 
uous though the uniformity is ensured. The width of 
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FIG. 12: Data collapse with a— 0.5, v = 1 



ject of our further research. To study mixed RFBM with 
WD, a quasistatic approach is employed to determine the 
critical stress which decreases linearly with the fraction 
of class A fibers. Exponents obtained in both the above 
cases fall under mean field universality class. The most 
interesting feature of the mixed models is the tunability 
of the critical stress with the system parameters. 

The behaviour of critical stress in a one dimensional 
RFBM with one sided LLS and strictly UD is also stud- 
ied using an exact method and the complication arising 
due to the uniform nature of distribution is emphasized. 
We have also looked at the breakdown phenomena of 
RFBM on a random graph with coordination number 
3. The critical stress is obtained using the finite size 
scaling method. The mean field nature of transition is 
established in the limit N — > oo. 



the forbidden region tunes the critical stress of the sys- 
tem and thus exhibits a rich critical behaviour. At the 
same time, restrictions on the parameters of the distri- 
bution for which the analytical method holds good are 
also clearly pointed out. The burst avalanche distribu- 
tion of this type of discontinuous threshold distribution 
is expected to be very interesting and that is the sub- 
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